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ABSTRACT 


The extension of Van Leer’s preconditioning procedure to generalized finite-rate chem- 
istry is discussed. Application to viscous flow is begun with the proper preconditioning 
matrix for the one-dimensional Navier-Stokes equations. Eigenvalue stiffness is resolved 
and convergence-rate acceleration is demonstrated over the entire Mach-number range from 
nearly stagnant flow to hypersonic. Specific benefits are realized at the low and transonic 
flow speeds typical of complete propulsion-system simulations. The extended precondition- 
ing matrix necessarily accounts for both thermal and chemical non-equilibrium. Numerical 
analysis reveals the possible theoretical improvements from using a preconditioner for all 
Mach number regimes. Numerical results confirm the expectations from the numerical 
analysis. Representative test cases include flows with previously troublesome embedded 
high-condition-number axeas. 

Van Leer, Lee, and Roe recently developed an optimal, analytic preconditioning tech- 
nique to reduce eigenvalue stiffness over the full Mach-number range. By multiplying 
the flux-balance residual with the preconditioning matrix, the acoustic wave speeds are 
scaled so that all waves propagate at the same rate, an essential property to eliminate 
inherent eigenvalue stiffness. This session discusses a synthesis of the thermo-chemical 
non-equilibrium flux-splitting developed by Grossman and Cinnella and the characteris- 
tic wave preconditioning of Van Leer into a powerful tool for implicitly solving two and 
three-dimensional flows with generalized finite-rate chemistry. 

For finite-rate chemistry, the state vector of unknowns is variable in length. Therefore, 
the preconditioning matrix extended to generalized finite-rate chemistry must accommo- 
date a flexible system of moving waves. Fortunately, no new kind of wave appears in the 
system. The only existing waves are entropy and vorticity waves, which move with the 
fluid, and acoustic waves, which propagate in Mach-number dependent directions. The 
non-equilibrium vibrational energies and species densities in the unknown state vector act 
strictly as convective waves. The essential concept for extending the preconditioning to 
generalized chemistry models is determining the differential variables which symmetrize 
the flux Jacobians. The extension is then straight-forward. 

This algorithm research effort will be released in a future version of the production- 
level computational code coined the General Aerodynamic Simulation Program (GASP), 
developed by Walters, Slack and McGrory. 
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ompressible Euler equations decouple at low speed. 
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